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ABSTRACT 

A ten-level primitive equation ocean circulation model 
is used to investigate the formation and. propagation of baro- 
clinic wavelike responses to annual variations in wind stress 
forcing. At 12.2°N a thermal response to Ekman pumping is 
observed to be in phase across the entire longitudinal extent 
of the model Pacific Ocean. A zonal fluctuation in interior 
ocean current systems is also suggested. Eastern boundary 
disturbances not related to local forcing are observed to 
propagate poleward. At 6.1°N, main thermocline displacements 
are observed to propagate westward as the result of vertical 
advection induced by Ekman pumping. Surface layer tempera- 
ture disturbances also exist which are produced by meridional 
advection of mean temperature (colder to the south) by sur- 
face Ekman currents. Eastern boundary disturbances also 
exist at this latitude and extend to a deeper level. All 
disturbances observed in this investigation are considered 
to be permanent features of the model circulation. Grid size 
limitations preclude detailed investigation of the observed 
boundary disturbances with this model . 
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INTRODUCTION 



I . 

Recent modeling investigations of the strength and 
thermal structure of Tropical North Pacific Ocean currents 
have provided evidence as to the existence of a baroclinic 
response to annual variations in wind stress. Meyers (1978) 
has observed variations in the depth of the 14°C isotherm 
which, at 6°N, propagate westward at nearly the speed of 
free, non-dispersive Rossby waves. His model mechanisms at 
6°N are Ekman layer divergence (Ekman pumping) as a response 
to local wind stress, and planetary, geostrophic divergence 
inherent on a rotating sphere. At 10°N .Meyers observed 
therraocline depth variations which were nearly in phase 
across the Pacific with the Ekman pumping velocity calcu- 
lated from the local wind stress. White (1977) found an 
agreement in phase only at 150°W, and suggested a baroclinic 
response to annual forcing in the form of a long wave which 
results from the superposition of baroclinic Rossby waves 
emanating from the eastern boundary, and the local forced 
response. This resulting long wave propagated at twice the 
speed of free Rossby waves. 

The purpose of this study was to investigate further 
any low frequency responses to wind forcing in the North 
Equatorial Pacific and to consider whether any such responses 
are consistent with theory. The numerical model used is that 
of Haney (1974) with subsequent improvements involving 
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parameterization of convective mixing and surface wind stir- 
ring (Haney and Davies, 1976) and nonlinear eddy viscosity 
(Haney and Wright, 1975). This model is based on the hydro- 
static and Boussinesq approximations and is a 10-level prim- 
itive equation model in a closed rectangular basin with a 
uniform depth of 4 kilometers. The model basin extends from 
the equator to 65°N and westward for 90° of longitude from 
125°W. Horizontal resolution is approximately 3° of longi- 
tude and 2° of latitude with 33 uniformly spaced grid points 
in both the east-west and north-south directions. The ten 
vertical levels are placed at 10, 32.5, 62.5, 102.5, 162.5, 
462.5, 900, 1700, and 3100 meters below the undisturbed sea 
surface. The model has been developed and improved in order 
to study large-scale thermal anomalies observed in the 
Central North Pacific Ocean by the North Pacific Experiment 
(NORPAX) . 

The model was "spun up” for a number of years of 
simulated time sufficient to allow circulations and thermal 
patterns to establish themselves and then detailed analysis 
of the transient response in the tropical regions was made. 
The wind stress prescribed as a forcing function in the model 
from the equator to 30°N was that used by Meyers (1978) and 
described in detail by Wyrtki and Meyers (1975). North of 
30°N wind stress was calculated from geostrophic wind data 
obtained from an NCAR data tape of the Northern Hemisphere 
climatological atlases of Jenne e^ > (,1969a, b). The 
major region of interest was south of 20°N. This location 
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was where Meyers and White observed wavelike disturbances 
and it was reasonable to expect similar waves in this study. 
By analyzing any observed waves as to their phase relation- 
ship with the local wind stress curl, their phase speed, 
and their wavelength, possible mechanisms for their genera- 
tion and propagation may be suggested and compared with 
theory . 

By better understanding how the oceanic circulation 
changes as a result of variant wind forcing, more compre- 
hensive studies of large-scale air-sea interaction and 
climatic changes may be undertaken. 
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II. DATA SOURCES AND PROCESSING 



Data from two different sources, one covering the region 
north of 30*^N, and one covering the region south of SO'^N, 
were combined to produce a realistic wind stress field which 
was used as a forcing function in this modeling investiga- 
tion. The first data base was that processed and explained 
in detail by Wyrtki and .Meyers (1975). This monthly aver- 
aged wind stress data, sorted into 2° latitude by 10° longi- 
tude fields encompassing 29°S to 29°N and 125°E to 75°W, 
was generously provided by Dr. Steve Piascek of NORDA. Only 
the data which was appropriate for use by the model, i.e. 
that data from the equator northward and from 145°E to 125°W, 
were used in this study. Wyrtki and Meyers calculated these 
stress data for each. observation according to the formula 








V 



) 



where t is the wind stress, rectangular component positive 
to the east (X) and north (4)), p„ is the constant surface air 
density of 1.2 kg/m , c^ is the drag coefficient assumed con- 

_3 

stant in the trade wind zone at 1.5 x 10 , u„ and v„ are 

’ a a 

the surface wind components, and V = (u + v )^, the 

3 . 3 . 3 . 

magnitude of the surface wind. A quadratic interpolation 
scheme was used to transfer the Wyrtki-Meyers monthly wind 
stress component arrays, e.xtending from the equator to 29°N, 
to the model gridpoints. 
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The second source of data was an NCAR tape of monthly 
mean climatological surface geostrophic wind component data 
CJenne e;t , 1969a, b). These wind components had been 
sorted into 5° latitude by 5° longitude fields encompassing 
the model latitudes and ranging in longitude from 120°E to 
80°W. These data were used to calculate surface wind stress 
using the same formula as above. However, an atmospheric 

_3 

drag coefficient of 1.3 x 10 was used (Businger e;t al . . 
1971) and a 10% reduction in magnitude from geostrophic was 
assumed. The resulting stress components for the model 
domain north of 30*^N were interpolated to model gridpoints. 

By combining the stress component fields computed and/or 
interpolated from the two data sources, 33 x 33 fields en- 
compassing the model's geographical domain for individual 
months were obtained. The Wyrtki-Meyers stress data, in- 
cluded from the equator to 29°N, is based on actual wind 
observations in the tropics, while the NCAR data, from 30*^N 
to the northern boundary, is based on geostrophic wind cal- 
culations from monthly mean sea-level pressure data. It is 
felt that the resulting stress field represents the best 
total stress pattern that is presently available from monthly 
mean data over the North Pacific Ocean. 

Although the major feature of the two stress fields 
matched quite well at the common latitude of 30°N, a smooth- 
ing technique was employed to remove small scale irregular- 
ities and improve horizontal consistency. A weighted three- 
point meridional average given by 
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w. 



was selected with the weight, W . , a Gaussian function of the 

latitude (i is the longitudinal gridpoint index and j is the 

latitudinal index). W. is detailed and illustrated in 

3 

Figure 1. At the interface between the two data sources, 
30°N, W. approaches one, and at the latitudinal boundaries 

«J 

of the model it approaches zero. Essentially, then, at the 
interface the smoothing process was a i, J meridional 
filter, and near the boundaries was negligible. The merid- 
ional smoothing was performed twice with boundary values 
being unchanged. Additionally, a i, |, J monthly time 
filter was used once over the entire fields and is given by 



s 

T . 



i, j 



= .25rx . 



i, j ,k-l 



+ 2x. 



i, j ,k 



+ X . 



i, j ,k+l 



), 



where k indicates the month. 

The ocean circulation model (Haney e;t aji^. , 1978) was 
integrated over a total of nine complete annual cycles start- 
ing from the results of the previous model. The only change 
from the previous model is the introduction of the new wind 
stress fields as described above. After nine years of inte- 
gration it was clear that the upper layers of the ocean 
model, i.e. those most strongly affected by the annual 
variation in wind forcing, were essentially repeating them- 
selves one year after another. The results analyzed below 
are those from the eighth and ninth years of the nine year 
integration. Primary attention is devoted to the annual 
variation of model variables, less to the annual mean. 
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III. OBSERVATIONS AND RESULTS 



The trade wind field over the Pacific Ocean is one of 
the largest and most consistent wind fields on earth although 
little is actually known about its time and spacial fluctua- 
tions (Wyrtki and Meyers, 1976). Wind-driven ocean circula- 
tion studies are largely dependent on wind stress field 
fluctuations since the curl of the wind stress, and thus 
Ekman surface layer divergence or pumping, is a very crucial 
forcing parameter. Accordingly, maps of the surface wind 
stress and of its curl were analyzed for significant fluc- 
tuations which may be related to observable oceanic 
disturbances . 

Figure 2 is a depiction of the annual mean wind stress 
vectors over the model domain. The region of primary inter- 
est in this investigation is that south of 20'^N. As shown, 
the Intertropical Convergence Zone (ITCZ) has an annual 
mean position at about 8°N. Wyrtki and Meyers (1975) have 
shown that the curl of the wind stress and thus the Ekman 
pumping velocity is relatively insensitive to changes in wind 
strength and depend mostly on the locations of strong merid- 
ional shear, which show latitudinal variation throughout the 
year. In the northern hemisphere winter and spring, the 
northeast trades are at their strongest and the southeast 

trades at their weakest. The ITCZ lies about 3*^ south of 
its annual mean position and as seen in Figure 3, a map of 
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wind stress curl for a typical spring month, regions of 
strong positive wind stress curl appear between 5°N and lO^N. 
Two zones of particularly strong positive curl are apparent, 
one between 155°E and 175°E and another between 155°W and 
140°W. In the northern hemisphere summer and fall, the 
ITCZ shifts northward to about 13°N and the zone of strongest 
cyclonic curl is located between 10°N and 15°N as seen in 
Figure 4. 

The strongest meridional shear in the surface wind stress 
is observed in a band slightly north of the ITCZ which passes 
over the annual mean position of the ITCZ twice per year. 
Therefore the amplitude of the first annual variation of the 
shear is largest at the latitudinal extremes of the ITCZ and 
small at the annual mean position. Meyers (1978) computed 
the amplitude of the first annual harmonic of the Ekman pump- 
ing velocity and the depth of the 14°C isotherm in the North 
Pacific and found that both reached maxima at 6°N and at 10°N. 
Not surprisingly, these latitudes coincided with the zones of 
strongest cyclonic wind stress curl. As mentioned previously, 
at 10°N the Ekman pumping velocity and the depth of the 14°C 
isotherm fluctuated in phase across the Pacific while at 6°N 
the temperature fluctuations propagated westward from a 
region of resonant forcing in the Eastern Pacific. In order 
to analyze results which could best be compared with those 
of Meyers (1978), model latitudes of 6 . 1°N and 12 . 2°N were 
chosen for study in this investigation. 
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a. Results at 12.2°N 

Meyers' (1978) result at 10°N was an oceanic tempera- 
ture response which was nearly in phase (a two month zonal 
lag) with the local Ekman pumping velocity across the entire 
Pacific Ocean. White (1977) theorized that the annual fluc- 
tuation at any interior ocean location from 10°N to 20°N 
depends on the fluctuating wind field over the entire east- 
ward ocean with the phase speed and amplitude strongly related 
to distance from the eastern boundary. In a bounded ocean, 
the presence of the eastern boundary can significantly alter 
the phase distribution of the oceanic response to Ekman pump- 
ing by acting as a concentrated source of baroclinic 
planetary waves which propagate boundary responses westward 
into the interior, drastically altering the phase relation 
predicted by local forcing only (White, 1977). 

.Model level 5, corresponding to a depth of 162.5 
meters, was chosen for detailed study since it is the model 
level which most closely corresponds with the observed depth 
of the 14°C isotherm (Meyers, 1978). This level is above the 
intermediate layers where the vertical temperature gradient 
is weak and below the level of the maximum density gradient, 
whose stability shields it from temperature fluctuations 
associated with sea surface heat exchanges. Time-longitude 
maps of departures from the annual mean (denoted with a prime) 
of the zonal (u') and meridional (v') flow, vertical velocity 
(w*), and temperature (T' ) were developed to show any zonal 
propagation of wavelike disturbances. Figures 5 and 6 are 
such maps of u' and v' respectively at level 5 across 12.2°N 
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over the last two years of the model run. These figures in- 
dicate westward propagating annual signals which emanate from 
near or at the eastern boundary with the zonal component hav- 
ing a slightly faster phase speed and a higher magnitude. 
Figure 7 is a similar depiction of the disturbance vertical 
velocity, w' . The major features of the open ocean vertical 
velocity field in Figure 7 develop near the eastern boundary 
and seem to propagate westward as far as the mid ocean region 
at which point further westward propagation with increasing 
time is not observed. White (1977) observed that due to the 
inverse latitudinal dependence in zonal phase speed, baro- 
clinic long waves in the eastern ocean are directed westward 
while in the western ocean they appear to be directed more 
northward due to refraction effects. Regardless, the overall 
pattern observed in Figure 7 does not correlate with those 
of the zonal and meridional flow components shown in Figures 
5 and 6 respectively. Vertical velocity is generally related 
to the divergent part of the flow; the propagating features 
of Figures 5 and 6 are most likely due to the rotational part 
of flow and as such have little impact on the vertical 
velocity. It is noted that disturbances in Figures 5 through 
7 appear to intensify along the western boundary of the model 
ocean. This is to be expected since westward propagating 
features tend to congregate at or near western boundaries 
in models such as that used in this study. 

Figure 7 also shows local disturbances of great inten- 
sity along the eastern boundary. However, these disturbances 
do not appear to propagate westward to any significant extent. 
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They also do not appear to be related to the longshore (merid- 
ional) component of wind stress nor to local Ekman pumping. 

In order to investigate any meridional propagation of these 
boundary features, a time-latitude map of w' along the eastern 
boundary was prepared and is shown in Figure 8. Wyrtki (1975) 
has surmised that the annual baroclinic response at an eastern 
boundary away from the equator is in part composed of a per- 
turbation which propagates into the region from the equatorial 
region in the form of an internal Kelvin wave. Those coastal 
disturbances at 12.2°N in Figure 8 have a computed (northward) 
phase speed of 11.7 cm/s and the major features have an annual 
period. It is interesting to note that Enfield and Allen 
(1978) have recently observed spatially coherent disturbances 
in monthly mean anomalies of temperature and sea level along 
the Pacific coast of North and South America. In the tropical 
latitudes, these observed fluctuations are also uncorrelated 
with local forcing mechanisms and propagate poleward at phase 
speeds comparable to those obtained in the model (Allen, 
personal communication). Because of the coarse grid in the 
present model however, no further comparison between model 
results and coastal observations is justified. 

Another important re.sult which can be seen from Figures 
5 through 8, and others below, is that the two consecutive 
model years are virtually identical. This means that the 
model is in a statistical equilibrium with the forcing and 
that the oscillations being examined are not transitory but 
are permanent features of the model circulation. 
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Figure 9 is a time-longitude map of temperature dis- 
turbances which is very consistent with the vertical velocity 
pattern given in Figure 7. A two month phase lag is apparent 
with generally warm water sinking during the first half of 
the year and cold water rising during the second half - 
clearly and indirect, thermally forced disturbance. Calcula- 
tions based on the annual positive disturbance in the east- 
central ocean in Figure 9 revealed a zonal wavelength of 
about 10,000 kilometers and a zonal phase speed of approxi- 
mately -42 cm/s. These results are in qualitative agreement 
with White (1977) who theorized a wavelength of 10,000 kilo- 
meters and a phase speed of -32 cm/s and observed a wavelength 
of 16,000 kilometers and a phase speed of -46 cm/s at 12^N. 

The Ekman pumping velocity was calculated from the 
following formula as described in detail by DeWitt and 
Leetmaa (1978) 



0 



W 



-h 
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A map of the disturbance Ekman pumping velocity is presented 
in Figure 10. This pattern corresponds with the gross fea- 
tures of the vertical velocity field of Figure 7 although no 
propagation is noted in the eastern ocean. This is in con- 
trast with the results of both White (1977) and DeWitt and 
Leetmaa (1978) who noticed a tendency for the Ekman pumping 
disturbances to propagate westward. Magnitude differences 
between corresponding features in Figures 7 and 10 are perhaps 
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due to the fact that Figure 10 represents the depth of the 
bottom of the Ekman layer which need not be coincident with 
level 5 of the model. Additionally, different processes 
are involved in vertical velocity and Ekman pumping. The 
theoretical Ekman depth was calculated according to the 
fonnula 
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(DeWitt and Leetmaa, 1978) where is the vertical eddy vis- 

2 

cosity coefficient in the model which is 10 cm /s, and f is 
the Coriolis parameter. At 12.2°N, this depth was calculated 
to be approximately 20 meters. 

In order to describe the vertical structure of the 
annual wave disturbances, vertical cross sections (_x,z) of 
the disturbance quantities were prepared across 12.2°N on 
August 3 of the last model year. This date was selected so 
as to coincide with the large positive amplitude disturbance 
noted on the time-longitude section of Ekman pumping in 
Figure 10. Figure 11, the vertical cross section of u on 
August 3, may be compared with the time-longitude map of u 
in Figure 5. The positive wavelike disturbance near the 
eastern boundary in Figure 5 extends from the surface to the 
bottom with greatest magnitude near the surface, as is the 
case with the positive disturbance in the western ocean. The 
negative disturbance in the mid-ocean extends from the surface 
to about 1500 meters with greatest magnitude at about 50 meters. 
The result that the largest magnitudes are not at the surface 
in Figure 11 may be explained by the dependence of the zonal 
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flow on temperature and pressure gradient structure and not 
so much on the surface winds. Generally, then, the zonal 
disturbances extend quite deep at their formative stage and 
do not significantly deepen or shoal as they propagate west- 
ward . 

Figure 12 shows the corresponding cross section of 
the meridional velocity perturbation which may be compared 
with the time-longitude section of Figure 6. The positive 
disturbance in the mid-western ocean in Figure 6 on August 3 
extends vertically to within about 50 meters of the surface. 
The negative disturbance in the eastern ocean extends from 
about 1500 meters to the surface where it expands to cover 
the entire model domain across 12.2°N. The negative distur- 
bance is of greatest magnitude near the surface at all longi- 
tudes while the positive disturbance is very small whenever 
it appears. Apparently the meridional flow component is 
more closely related to the surface winds than is the zonal 
component . 

Figure 13 depicts the vertical-zonal cross section 
of vertical velocity disturbance across 12.2°N on August 3. 
The positive disturbance in Figure 7 encompasses most of the 
model domain except for the surface layer and the extreme 
eastern boundary. Largest magnitudes are centered at depths 
of about 1000 meters. The eastern boundary negative distur- 
bance extends from about 200 meters to the surface. 

Corresponding temperature disturbances are shown in 
Figure 14. Warm conditions are noted in all areas on August 
3 except the eastern region at levels below about 75 meters. 
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Wunsch (1977) has noted that if wind forcing is a.s- 
sumed to be over an annual period then there is no reason to 
isolate the problem of vertical propagation from that of 
horizontal. He contends that internal gravity waves will 
determine the vertical modal structure and that in the 
presence of a non-dissipat ive and flat bottom energy may 
travel from the surface to the bottom and then be reflected 
back to the surface. Any vertical propagation of wavelike 
disturbances should be apparent in time-depth maps of the 
perturbation quantities. These depictions were developed 
and are shown in Figures 15 through 18. These maps were made 
at an arbitrarily selected interior ocean longitude of 162°W. 
No vertical propagation is seen in the zonal and meridional 
velocity components, Figures 15 and 16 respectively. Again 
it can be seen that u' penetrates deeper than v' and is 
therefore most likely due to the north-south meandering of 
the major equatorial current systems which are primarily 
zonal. Figure 17 reveals that the vertical velocity distur- 
bances also do not appear to propagate vertically but they 
attain their greatest magnitudes at fairly substantial 
depths. These deep centers apparently are the cause of the 
corresponding deep temperature disturbances (Figure 18) 
which lag the vertical velocity by about 2§ months. The 
temperature fluctuations above 150 meters are due to the 
seasonal cycle of surface heat flux modified slightly by a 
fluctuating meridional temperature advection in the upper 50 
meters. Again notice how the two model years are nearly 
identical - thus in equilibrium with the forcing. 
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The model results at 12 . 2°N are summarized as follows 
The dominant thermal response at the depth of the main thermo 
dine occurs in phase across the entire longitudinal extent 
of the ocean and is attributed to large scale vertical advec- 
tion produced by wind-induced Ekman pumping. This is in 
agreement with the theoretical analysis and observational 
results of Meyers (1978). In addition, there is a pronounced 
fluctuation in the zonal current in the open ocean which 
appears to be associated with a north-south seasonal meander- 
ing of the main zonal current systems. There also exists 
poleward propagating disturbances confined to the immediate 
vicinity of the eastern boundary which do not appear to be 
related to the local forcing. Further study is required to 
determine whether these disturbances are in any way related 
to the sea level disturbances observed along the North 
American coast by Enfield and Allen (,1978). 
b. Results at 6 . 1°N 

o 

Variations in the depth of the 14 C isotherm which 
propagated westward at nearly the phase speed of free non- 
dispersive Rossby waves were observed at 6°N by Meyers (1978) 
He theorized that the waves were generated by the zonal and 
seasonal variability in Ekman layer divergence in the 
eastern Pacific and that they propagated nearly freely in 
the western Pacific where he found the forcing was weak. The 
resulting wave lagged the Ekman pumping velocity by slightly 
more than one month at the position of maximum forcing. 
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Figures 19 and 20 are time-longitude maps of zonal 
and meridional flow disturbances, u^ and v' respectively, 
across 6.1^N at 162.5 m (model level 5) during the same two 
model years analyzed above at 12.2°N. These figures indicate 
rapidly westward propagating annual signals which appear to 
emanate from the eastern boundary. The zonal component in 
Figure 19 appears to increase gradually in magnitude to a 
maximum in the region slightly west of the mid-ocean and then 
decreases from there to the western boundary. The zonal dis- 
turbances propagate faster than the meridional and faster 
at 6 . 1°N than at 12 . 2°N (Figure 5). Figure 21 is a time- 
longitude map of vertical velocity disturbance at the same 
latitude and depth. As was noted for Figure 7 at 12.2°N, 
local disturbances of great magnitude are observed along the 
eastern boundary which do not appear to propagate westward. 
These northerly propagating Kelvin-like disturbances (it can 
be seen from Figures 19 through 22 that v' , T' , and w' are 
involved in the disturbance but not u') are depicted in 
Figure 8. Near the eastern boundary in Figure 21, wavelike 
perturbations form and propagate completely across the model 
domain, in contrast with those at 12.2°N (Figure 7) which 
appear to be nearly in phase across the ocean. The distur- 
bances' amplitudes slowly vary with westward propagation. 

The temperature disturbance at 6 . 1°N (Figure 22) 
has essentially the same pattern as the vertical velocity 
disturbance but with a two month phase lag which indicates 
the temperature is responding to the vertical velocity. Cal- 
culations revealed a phase speed of -32 cra/s, a zonal 
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wavelength of about 9900 kilometers, and an annual period. 
Meyers (1978) observed a corresponding phase speed of -56 
cm/s with a wavelength of 17,000 kilometers at 6°N. 

A time-longitude map of Ekman pumping velocity distur- 
bance across 6 . 1°N is shown in Figure 23. The theoretical 
Ekman depth as calculated above for 12.2°N was found to be 
approximately 28 meters at 6.1°N. The pattern in Figure 23 
correlates well with the actual vertical velocity (Figure 21) 
which in turn correlates with the temperature fields except 
for the region of very slow westward propagation in the mid- 
ocean area. Propagation is very rapid in both the eastern 
and western thirds with greatest amplitudes observed at the 
eastern boundary (125°W). The latter observation is consis- 
tent with Meyers' (1978) finding that the maximum displace- 
ment occurs at 130°W. The region of strong cyclonic wind 
stress curl observed at this latitude in Figure 3 was at 
and slightly west of 140°W. 

The zonal-depth structure of the disturbance patterns 
across 6.1°N on August 3 of the last model year is depicted 
in Figures 24 through 27. Figure 24 shows that the positive 
u' disturbance of Figure 19 is very broad at great depths, 
narrow at mid levels, and widens again in the surface layers 
with slightly higher magnitudes observed at the surface. The 
corresponding negative disturbance is widest and has its 
greatest magnitude at about 75 meters. As was discussed pre- 
viously with regard to Figure 11, there is a deep zonal 
current fluctuation which seems to be related to the north- 
south displacement of the main current systems which are 
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zonal . 



The meridional velocity disturbance across 6 . 1°N is 
presented in Figure 25. A relatively weak flow is predomi- 
nant across the model domain except in the surface layers 
where, as in Figure 12, there is a strong response due to 
the dependence on surface winds. Positive disturbances shown 
in Figure 25 are of much smaller magnitude. 

The zonal-depth map of w' in Figure 26 shows that 
with the exception of a narrow semi-vertical band in the 
western ocean and the Kelvin-like disturbances on the eastern 
boundary, the model domain is predominantly occupied by sink- 
ing (w' < 0) disturbances at this time. The shallow Kelvin 
disturbance extends deeper at 6.1°N than it did at 12.2°N 
(Figure 13) and a phase change is apparent between the two 
latitudes. Figure 27 shows that the model domain is colder 
than normal except for the eastern ocean surface layer and 
a mid-level warm disturbance in the western ocean. 

Time-depth maps of u' and v' at 162°W are shown in 
Figures 28 and 29. Disturbances of greatest magnitudes occur 
at the surface with those of the zonal component e.xtending 
quite deep. On comparing these figures with their counter- 
parts at 12.2^N, Figures 15 and 16, the differences in dis- 
turbance magnitudes between the two latitudes is obvious; the 
deep response is a little stronger at 6 . 1°N but the surface 
Ekman flow is nearly an order of magnitude higher at 6.1°N. 
Plots of the disturbances of the zonal and meridional com- 
ponents of wind stress as functions of time at 162°W are 
shown in Figure 30. The large magnitude surface layer distur- 
bances noted in Figures 28 and 29 are in phase with these 
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stress disturbances indicating that the wind stress is most 
likely the dominant forcing process of ocean (Ekman) flow 
at the surface and immediately below. The time-depth map 
of vertical velocity disturbance across 6 . 1°N as presented in 
Figure 31 is similar to that of 12.2*^N (Figure 17) in that 

greatest magnitudes are reached at deep levels. However, in 
contrast to 12.2°N, there is a phase shift at intermediate 
levels at 6.1°N. The phase shift probably represents a 
transition region between the upper 500 meters or so, which 
is responding to the wind, and the remainder of the ocean 
which is reacting out of phase. Figure 32 shows the strong 
response of the temperature field in the upper 500 meters to 
the vertical motion with a two month lag. At the same time, 
surface disturbances due to seasonal wind changes occur but 
they are confined to the upper 50 meters. At 6 . 1°N v' advects 

cold water northward from the equatorial surface region during 
the first half of the year and warm water southward during the 
second half of the year. 

The model results at 6 . 1°N are summarized as follows. 
The dominant thermal response in the vicinity of the main 
thermocline (level 5, or 162.5 meters in the model) is a west- 
ward propagating disturbance (Figure 22) produced by vertical 
advection (w' in Figure 21) induced by Ekman pumping (Figure 
23). Figure 27 shows that the temperature response seen at 
162.5 meters involves the entire thermocline and is actually 
strongest at about 100 meters. In addition to the propagating 
thermal response in the main thermocline there is an equally 
pronounced annual cycle of surface temperature Cl'igure 32) 
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which is produced primarily by meridional advection of mean 
temperature (colder to the south) by surface Ekman currents 
(Figure 29). 
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IV. CONCLUSIONS 



Wavelike features simulated in the tropical regions of 
an ocean circulation model which are seen as disturbances 
from the annual mean of vertical velocity and temperature 
have been described. These features are considered to be 
permanent oscillations of the model circulation and result 
from variability in wind stress forcing. At 12.2°N, a 
thermal response to large scale vertical advection produced 
by wind-induced Ekman pumping or surface layer divergence 
occurs practically in phase across the entire longitudinal 
extent of the model ocean. Underlying features of the 
annual signal do appear to propagate westward. These results 
are in agreement with the theory and observations of .Meyers 
(1978). Evidence also suggests a zonal fluctuation in the 
interior ocean current systems. Along the eastern boundary 
poleward propagating perturbations of the vertical velocity 
field exist which may be related to sea level disturbances 
which have been observed by Enfield and Allen (1978). 

At 6 . 1 '^N, thermocline displacements propagate westward 
also as a result of vertical advection induced by Ekman pump- 
ing. Disturbances at this latitude are significantly higher 
in magnitude than those observed at 12.2°N. Additionally, 
pronounced sea surface temperature disturbances are apparently 
produced by meridional advection of mean temperature by 
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surface Ekman currents. Eastern boundary perturbations also 
exist and extend deeper at 6.1°N than at the higher latitude. 

Results and observations in this study were generally 
consistent with accepted theory. It would be beneficial to 
consider latitudes other than 6 . 1*^N and 12.2^^N in future 
studies in order to more completely describe the tropical 
wavelike disturbances which were noted. The equatorial 
boundary in the current model may, however, cause unreason- 
able results at lower latitudes than 6.1°N. The model is 
presently not suitable for use for further study of the 
eastern boundary disturbances because of grid size limita- 
tions. However, a careful study of these features using a 
higher resolution model having a realistic coastline con- 
figuration may be valuable in light of the recent findings 
of Enfield and Allen (,1978). 

The present model has little or no synoptic ’’noise". 

It is recognized that this feature may be altered if finer 
resolution or synoptic wind forcing were introduced. 
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